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Abstract 

A special case of the Fateev-Zamolodchikov model is studied resulting in a solution of the Yang- 
Baxter equation with two spectral parameters. Integrable models from this solution are shown 
' to have the symmetry of the Drinfeld double of a dihedral group. Viewing this solution as a 

descendant of the zero-field six-vertex model allows for the construction of functional relations 
and Bethe ansatz equations. 



^ ■ 1 Introduction 
>■ . 

. The connetion between the Yang-Baxter equation and integrable systems has been well-studied 

' with origins in the works of McGuire, Yang and Baxter [5, 31, 41]. Solutions of the Yang-Baxter 

equation, i?-niatrices, allow the construction of integrable one-dimensional quantum chains via the 
quantum inverse scattering method (QISM) [40]. Typically the quantum chains produced consist of 
C^l . nearest neighbour interactions, are translationally invariant and have periodic boundary conditions. 

However, it is possible to modify this procedure resulting in a model with twisted boundary conditions 
[14, 21, 37], in which periodicity is retained but translational invariance is broken as the interaction 
between the first and last site is distinct from interactions between other neighbours. The periodicity 
of the model can also be broken using Skylanin's approach [39]; this procedure was termed BQISM 
' (boundary quantum inverse scatter method) and leads to quantum chains with open boundary 

^ . conditions. Skylanin's approach required the advent of reflection matrices and assumed multiple 

properties of the i?-matrix involved. Over time these assumed properties have been systematically 
removed [9, 30, 32, 33, 42]. Operators which lead to a quantum chain with open boundaries can be 
modified to produce a quantum chain with braided closed boundary conditions [20, 24, 26, 28, 29, 27]. 

Connected to the Yang-Baxter equation are the algebraic structures known as quasi-triangular 
Hopf algebras. These structures are known to provide algebraic solutions of the Yang-Baxter equa- 
tion, albeit without so-called spectral parameters. The addition of spectral parameters can be 
achieved through Baxterisation [11, 25]. The Drinfeld doubles of finite group algebras are one no- 
table class of quasi-triangular Hopf algebras [15]. Although attention has recently been paid to these 
algebras in relation to non-abelian anyonic theories, not much interest has been directed towards the 
connection between their associated i?-matrices and integrable quantum chains. We remark that in 
contrast to the more familiar case that the i?-matrix comes from a quantum group, the quantum 
chain is interpreted as an anyon chain [10, 17] rather than a spin chain. 

We are concerned with the Fateev-Zamolodchikov model [16] , which is known to be a limiting case 
of the chiral Potts model [3] . The Fateev-Zamolodchikov model has received attention for both open 
boundary conditions [43] and periodic boundary conditions [1, 2, 36]. The latter works are restricted 
to the uniform square lattice case of the checkerboard lattice model [6, 35]. This special limiting case 
is found by making certain rapidities equal. We however are concerned with an alternate limit of 
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the Fatccv Zamolodchikov model. In this hniit we perform a basis transformation on the i?-matrix, 
again constructed as checkerboard vertex model, where the symmetry of D{Dn) is recognisable. We 
discuss three different classes of boundary conditions for quantum chains, namely, periodic, open 
and braided closed boundary conditions. Lastly wc present the i?-matrix as a descendant of the 
six-vertex model and construct some fusion relations for the models discussed. 



2 Preliminaries 

Throughout §1 will represent the Kronecker delta fimction, being one if i. = j and zero otherwise. 
We do not restrict i and j to the integers and they may belong to any group. 

We define Cij to be an elementary matrix of Mdxd(C) whose indices are considered modulo d. 
These matrices obey the relation 



j = k mod d, 
j ^ k mod d. 



Our convention is that e^.j corresponds to the matrix with a one in the zth row and jth column and 
zeros elsewhere for 1 < i,j < d. For such operators we define the set of d-dimensional vectors {vk}k 
satisfying the property 

ej.kVk=Vj, Vj, fceZ. 

Under this definition we have that the indices of the vectors are considered modulo d as is the case 
with the elementary matrices. 

The Yang— Baxter equation 

The Yang-Baxter equation (YBE) is a non-linear equation in the three fold tensor product of some 
algebra A. There are multiple variants of this equation. The first we consider is the constant YBE, 
given by 

R12R13R23 = -^23^13-^12, 

where the indices of the operator R G A^A indicate the spaces which the operator acts non-trivially 
on. That is, 

R12 = R(g) I, R23 = I (g) R, etc., 
where / is the identity element of A. The parameter dependant form of the YBE is given by 

Ri2{x;y)Ri3{x; z)R23{y; z) = R2z{y]z)Riz{x;z)Ri2{x\y), (1) 

where .t, y and z are elements of the Cartesian product of d copies of C, which we consider a group 
under multiplication. This equation can be simplified if the i?-matrix satisfies the multiplicative 
analogue of the so-called difference property, 

Ri2{x)Ri3{xy)R2M = R23{y)Ri3ixy)Ri2ix). (2) 

Solutions to any of these equation will be referred to as R-matrices. 

When the i?-matrix is a matrix operator, i.e R{z) € End(y V) for some vector space V, we 
use the definitions of regularity and unitarity, 

R{i)=P and i?i2(5)i?2i(5"^) oc / (g) /, 

where P is the usual permutation operator 

P{v ^ w) = W ^ V, VjW gV. 

Wc assume every matrix operator i?-matrix satisfies regularity, unless otherwise stated. It is known 
that regularity implies unitarity. 
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The Drinfeld doubles of dihedral groups 

The algebraic structures which we focus upon are the Drinfeld doubles of dihedral groups. The 
general structure, along with the representation theory, of the Drinfeld doubles of finite group 
algebras is known [15, 23]. For any finite group, G, the defining relations of its Drinfeld double, 
D{G), can be presented purely in terms of the operation of the group. There are procedures which 
allow the construction of all irreducible representations (irrcps) of D{G). 

We arc interested in the dihedral groups of order 2n, where n is odd. The dihedral group i?„ 
corresponds to the symmetries of a regular polygon with n vertices. It has the presentation 

Dn = {f, T|cr" = = arar = e}, 

where e is the identity element of the group. We now consider the Drinfeld double of this group. It 
is defined as the vector space: 

D(D„)=C{gh*\g,h€Dr,}, 
equipped with the multiplication and coproduct 

giKg2h; = SlUll] {gig2)K and l^{gh*) = ^ gik-^h)* ® gk*. 

keD„ 

The remaining Hopf structure is uniquely determined by these relations. This algebra is associated 
with a canonical element, 

^= ^9^9*, 

geD„ 

which is known to satisfy the YBE. We define the twist map and twisted coproduct to be, 

T{a<^b) = b(^a and A'^(a) = T o A(a), \la,h & D{Dn), 
respectively. The ii-matrix satisfies 

9iA(a) = A^(a)$H 

for ah a e 

The explicit irreps of D{Dn) have already been constructed [13], however, we shall only use one 

of dimension n. Firstly we view D„ as a subgroup of Sm the group of permutations on n objects. 
This affords D„ an action on Z„, the ring of integers modulo n. This action is defined by 

n-l 
2 

CT = [1, 2, n] and r = JJ [A;,n — fc], 

fe=i 

using the standard notation of Sn [22]. Thus a increments an integer by one and r maps an integer 
to its negative, both modulo n. This action of £>„ allows us to present the n- dimensional irrep in 
the following manner: 

n n—1 

^(S') = H ^9{i),i and Tr{g*) = ^ Sf'^ejj. 

i=l j=0 

This representation acts on the space C" and can be naturally extended to (C")®' through the 
definition 

Trf (a)^^f (AW(a)), Va € 

where A*^') is defined recursively by the relations 

A^'^) = (A ® id®('=-2))A('=-i) and A^^) = A, 
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for A; > 3, with id being the identity map. 
The Fateev— Zamolodchikov model 

The models we consider are special cases of the well-known Fateev-Zamolodchikov model [16]. The 
Fateev-Zamolodchikov model is defined by weights 

W{z\Q) = l, W{z\l)^Y[ X^3-i-z ^ 



W{z\0) = 1 and W{z\l) = J| 



. , X'iz-l ' 

where A is a primitive 2nth root of unity and z & C and < Z < n — 1. These weights satisfy the 
star-triangle relation^ and are extended by the relations 

W{z\l) = W{z\n + l) = W{z\-l) and W{z\l) = W{z\n + l) = W{z\-l), 

for I e N. Using these weights we have the Fateev-Zamolodchikov iJ-matrix 

n 

R{x;y) = ^ W{xiy^'^\h - c)W{x2yi'^\h - d)W{x2y^'^\a - d)W{xiy^'^\a - c)ea,b® ec,d. 

a.b,c,d—l 

The weights here have been configured in a checkerboard lattice configuration. The uniform square 
lattice arises when the two different rapidities are set equal, xi = x^ and yi — y^ [6]. 



3 Integrable n-state models 

An observation of Bazhanov and Perk was that the D{D^) model presented in [13] was a special case 
of the Fateev-Zamolodchikov model [7]. Moreover it was implied that all the i?,-matric;es derived 
from D{Dn) as presented in [18] should also be special cases as was shown to be the case for odd n. 
Here we show that there is in fact a more general limit of the Fateev-Zamolodchikov i?-matrix that 
has the underlying symmetry of D(_D„). 

We take a limit of the Fateev-Zamolodchikov i?-matrix and define 



R{z) = lim lim R{x;y), 

X2,y2^0 Xl,3/l-)-C!0 



where zi = ^ and Z2 = Explicitly we have 

n 

a,6,c,d=l 

This has the properties 

i?(l,l) = P and PR{zi,Z2)P = R{z2'^,z^'^), 

This is now an i?-matrix which satisfies the variant of the YBE given by Equation (2). To put the 
i?-matrix in a recognisable form we perform the basis transformation corresponding to the matrices 

^ n -. n 



See [35] for definition of the stax-triangle relation. 
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The ii-matrix becomes 



E 

a,i,j=l 



.6=1 



(3) 



where we have replaced A with —w ^, with w being a primitive n root of unity. We can express the 
i?-matrix in terms of the projection operators presented in [18]. The projection operators are 



I n-1 

-E 

i, (a,6) = (0,0), 

1, (a,6)^(0,0). 



where, 



The admissible pairs, (a, 6), for the projection operators along with the corresponding irrep they 
project onto can be found in Table 2 of [18]. Using these projection operators we can express the 
i?-matrix as 



PR{z) = £/(o,6)(5y°''') + EE/(«.^')(^)?'^"'' 



6=0 



a=l 6=0 



where 



/(a,6)(^) = 



,c=l 



.d=l 



(4) 



As the i?-matrix can be written as linear combination of the above projection operators in such a 
manner it follows that the ii-matrix intertwines the coproduct of D{Dn). We express this as 

[PR{i), 7r®^(a)] =0, Va e D{Dn). 

It is known that setting zi= Z2 = z recovers the ii-matrix presented in [18] . Looking at the functions 
.f{a,b) {z) we see that they factorise into functions of z\ and Z2 and that the components involving Z2 
arc determined purely by the difference a — 6. Both of these are required by the construction of the 
D{Dn) -R-matrix as a descendant presented in [18]. 
In this basis we also have the additional property 

[R{z^,Z2)]* =R{z;,zl). 

By investigation of the non-zero entries of the i?-matrix we sec that there is a generalised eight-vertex 
condition present. However, this generalisation differs from others [4] which stems from different 
conservation rules. The conservation rule which is satisfied here is 



Rf.(z)=Q when i + j^k + l (mod n) for the operator R{z) = ^ Rfj{z)ei^j ® ek,i. 

Important to constructing integrable quantum chains are the logarithmic derivatives^ of R{z). 
We define 



if (1) = 



d 

dzi 



MR{z)) 



and = -i 



J z=l 



dZ2 



HR{z)) 



2=1 



^We define the logarithmic derivative as: -j- ln[A{z)] = A~''-{z)-^[A{z) 
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Using these definitions we have 

n n — 1 



i,j=i 1=1 



i-iy 



w 



-2l{i-j) 



~2~ n — 1 



= ^EE 



a=l 6=0 



(-!)»- 



,2ab 



and, using ^ PH'^^)p^ 



n n—1 

^^EE 
ij=i i=i 

a=l 5=0 



(-ly 



„2K»-i) 



*EE(-ir 



-2a6 



E ^7(a-6),7(Ti-6) ® e^(a),7(n) 



w 



2a 2a 



6^(a-fc),7(n-b) 

7e-Dn 



These two operators are self-adjoint and complex conjugates 



^(1) ^ Tjd) 



and 



^(2) 



°7(a),7(") 



Wc have that the Eigenvalues of H^^^ and H^"^^ are real, furthermore, the Eigenspectrums of H^^^ 
and -ff must be the same. 

These operators will form the local Hamiltonians of the quantum chains discussed in the next 
section. Although it is possible to present them in terms of spin operators, there is no physical 
motivation for this as the models will lack the usual underlying spin symmetries. Instead we see 
the emergence of the D{Dn) symmetry in these operators with summation over the elements of Z)„, 
implying the known invariance. The coefficients of the fusions channels, indicating energetically 
favoured paths, can be calculated from the functions given in Equation (4) [34, 38]. 



3.1 Quantum chains 

Here we will discuss different models which can be constructed. To describe the models we consider 
a general framework, in which we have combined quantum chains with periodic, twisted, open and 
braided closed boundary conditions. For this we replace the single _R-matrix with four operators, 
noting that for each class of model considered the different operators will be trivially related to the 
original i?-matrix or themselves trivial. The four operators we consider are those appearing in the 
equations, which they must satisfy, below. 



Ri2{x)Ri3{xy)R23{y) = R23{y)Ri3{xy)Ri2{x), 

R[2{x)Ri3{xy)R23{y) = R23{y)Ri3{xy)R[2{i), 

Ri3{x)Ri2{xy)R23{y) = R23{y)Ri2{xy)Ri3{x), 

Rl3{x)R\2{xy)R23{y) = R23{y)R'i2{xy)Ri3{x). 



(5) 

(6) 

(7) 
(8) 



The first of these equations is the usual aforementioned Yang-Baxter equation while the last three 
are simple variants. We also introduce two reflection equations along with reflection matrices 

Ri2i£y-')K^{S:)Ri2{S;y)K^{y) - K^iy)R[2ixy)K^ix)R'i2ixy-'), (9) 

[R[,{xy-')r'K+{5:)R[^{S:y)K+{y) = K+ {y)R^2{iy)K+ {i) [R^2{5:y-')r\ (10) 



where 



Ri2{S) [Ri2iS)] =1(31 and R'-^^iz) [R[^iz)] ' = 1(3 1 . 



t2 



6 



Wc have assumed the invertibihty of i?*2(^) ''■nd which is not always guaranteed. Provided 

that the above six equations arc satisfied then it follows that the monodromy matrix 

To^-L..c){z) = Roc{~z)...Roi{z)K^{z)Roi{z)...Roc{~z), 

satisfies 

i?i2(5?}-i)Ti3(5;)i?i2(iy)r23(y) = r23(y)i?'i2(iy)T^i3(i)^'i2(iy"')- 



This leads the transfer matrix 



which is guaranteed to satisfy 



[t{i),t{y)] = ^- 

This commutation relation ensures the integrability of the Hamiltonian: 



dZ; 



Mm 



[dim(y)]"^ tr 



dZ; 



Mt{z)) 



where aj G M are free and 13 j G 



are fixed. The later scalars are chosen and fixed to ensure self- 
adjointness when possible. We have also included a term to make the trace of the global Hamiltonian 
zero. 

For each model we will be able to write the global Hamiltonians in terms of other Hamiltonians 
acting on the entire space. That is, 

d 

with each of the 'H*^-'' and 'H^'^) commuting. Instead of calculating the Eigenspectrum of "H, often 
the individual spectra of the H^^^ will be calculated. However, to recover the spectrum of the com- 
plete Hamiltonian so-called pairing or grouping rules will need to be found. Generally the ground 
state (and its energy) will depend upon the aj with level crossings and quantum phase transitions 
expected. 



Periodic and twisted boundary models 

To construct a periodic model we set 

R'{z) = R{z) and R{z) = R'{z) = / O 7. 

We immediately find Equations (6-8) are trivially satisfied. Furthermore without loss of generality 
we can set K~^{z) = I satisfying Equation (10). The only operator left to discuss is K~{z) which 
must satisfy 

R{xy-'){K'i'i) ® K-{y)) = {R- {i) ® R- {y))R{iy-^). 

This equation can clearly always be satisfied although the spectrum of solutions depends heavily on 
R{z). The transfer matrix simplifies down to 

t{z) = tro [K^{z)Roc{z)...Roi{z)\ . 

In the case that K~{z) is independent of z the model presented is precisely the definition of a 
periodic model with twisted boundaries as defined in [21]. For the i?- matrices we are concerned with 
there is a class of such solutions. For each group element g S we can set 

K-{z) = TT{g). 

The global Hamiltonian will be 

c-i 

■H = Y,Hiii+i)+H^ci (11) 
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where 

d 

i=i 

is the local Hamiltonian and 



■j-HR{z)) 



dz 



-tr[dim(y)]"^tr 



= {nig)-'^I)H{Tr{g)^I), 



is the twisted local Hamiltonian. Explicitly we have for the iJ-matrix given by Equation (3) that 
the local Hamiltonian is a linear combination of the operators H^^^ and H^'^\ 



H 



2 n — 1 

EE 



(-1) 



o=l 6=0 

while the twisted local Hamiltonian is 

n—1 



la 



-2a 



^l{a-b),^{n-b) 



■'7(a),7(n)i 



(12) 



^^ = ^EE 



a=l 6=0 



(-1) 



X] e^(a-6),7(n-6) 



°907(o)>9°7(")' 



The global Hamiltonian is self-adjoint. 

Setting g = e recovers a periodic chain, while additionally setting ai = 1 = —a2 yields the 
Hamiltonian for the i?-matrix described in [18]. Wc have the property that the local untwisted 
Hamiltonian commutes with both the coproduct and twisted coproduct, that is, 



[Hi2, {it (g) 7r)A(a)] = [H21, (tt (g) 7r)A(o)] = 0, 



Va e D{D„). 



However, wc find that for periodic chains of finite length the global Hamiltonian does not inherit the 
full symmetry of D[Dn)- In this case we find that the inherited symmetry of the global Hamiltonian 
is instead given by 

C7„ = {a e £>(£>„)! A(a)=A^ (a)}. 

That is, 

[H,(7rr)(«)]=0, aeC«. 
The periodic chain of any length will be translationally invariant as 

t{l) = P1C-P13P12, 

must commute with the Hamiltonian. 

For general g we can see that the inherited symmetry of the system is reduced further. The 
global Hamiltonian will not be translationally invariant. However, we can consider 

t{l)=K^{l)Pl2P23:PiC-l)C 

to play the role of a generalised translation operator. To now leave the model unchanged we need 
to perform a translation, i.e. map site i to site i + 1 modulo n, and additionally perform a trans- 
formation on the 1st site. We also remark that if we have group elements g and h which are in the 
same conjugacy class then the global Hamiltonians produced from both are equivalent. 

Open boundary models 

To construct an open model we set 

R[^{i) = Ri2{S) = R2i{S) and R'^^iS) ^ Ri2iS). 

We find that with these definitions Equations (5-8) are all equivalent and both (9,10) are non-trivial. 
One particular case is when the ii'-matrices are the identity i.e. 



K-{i) = K+{z) = I. 
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This is always possible for the ii-matrix presented. This gives the global Hamiltonian 

c-i 

i=l 

where H proportional to the local Hamiltonian presented in the periodic case. This system is 
described as having non-interacting boundary terms or free ends due to the lack of operators acting 
solely on sites 1 or C. Unlike the periodic case the global Hamiltonian inherits the complete symmetry 
of D{Dn) for finite size chains. The introduction of interacting boimclary terms will break this 
symmetry. We find the model has no form of translational invariancc and that the operator t{l), 
which we previously identified with translation, is a scalar multiple of the identity. The latter 
property of the model also implies that the logarithmic derivative used for constructing the global 
Hamiltonian becomes the usual derivative with a scaling factor, as is the case with standard BQISM. 

Other boundary conditions have been studied for the Fateev Zamolodchikov model. Boundaries 
were considered in [43] using the language of weights and the star-triangle relation. In the particular 
case of the original D{D^) limit iiT-matrices were studied in [12]. 

Braided closed boundary models 

To construct a braided model we set 

R[^(z) = R2i{z), R= Mm R'{z) and R' = Mm R{z), 

Z^Zo Z—^Zq 

where 

zo G {(0,0), (0,00), (00,0), (00,00)}. 

Using these we have Equation (5) implies Equations (6-8). Reflection matrices K^{z) and K~{z) 
are still required, however, given K^{z) and K^{z) which lead to an open model we can set 

= lim K^{z), 

Z^Zo 

obtaining an integrable model. Here we will consider the case where = I, this will lead to the 
global Hamiltonian 

c-i 

i=l 

where H is the local Hamiltonian described by Equation (12), 

G = t{i) = 6i62...6£-itro [PocRoc] and bi = Ri(i+i)Pi{i+i). 

For all choices of zq we have that tro [fozi-Roz;] is a scalar multiple of the identity and the global 
Hamiltonian is invariant under the action of _D(Z?„). Like the case with the; twisted boundary model 
we have that t{l) again plays the role of a generalised translation operator. Here t{l) is a product 
of solutions to the braid equation, 

hbi+ibi = bi+ibibi+i. 

Thus the model is invariant under translation through braiding. 
We have four possibilities for zq but from unitarity we know that 

Urn R'{z)P = hm [i?'(5-i)P] , 

z — yzo z — >zo 

thus the braiding operator 6j generated from (0, 0) will be the inverse of the braiding operator from 
(00,00). Likewise for the two remaining possibilities. If Zo = (0,0) then we find that the braiding 
operator comes from the representation of the canonical element associated with £>(D„), 

n 

bi = Pi2(7r 7r)fHi2 = ^ Ci^i+j (g) ei_j,j. 
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In the case that zq = (0, oo) the braiding operator is not of such an elegant form, 

n 

a,z,j,fe=l 

For some other models it is possible to set 61 = P, in which case periodic chains appear. 

We remark that in contrast to earlier work on braided models [28] we have not assumed any 
properties of the i?-matrix except for regularity. This modification was possible because the braided 
models can be viewed as modified open models and subsequently the original approaches retained 
many of Skylanin's unnecessary assumptions. If we further take this view then it appears regularity 
is also not necessary [42]. 



4 Descendant of the Six- Vertex model and Constructing Fu- 
sion Relations 

The general _R-matrix presented in [18] was constructed as a descendant of the zero-field six-vertex 
model while requiring the symmetry of D{Dn). Here we started with the Fateev-Zamolodchikov 
model, a known descendant of the zero- field six- vertex model, and took a special limit in which _D(Z)„) 
symmetry emerges. We use this connection to the six- vertex model to derive fusion relations. 
The zero-field six-vertex ii-matrix is given by 



r{z) 











V 



w 



w 



\ 



^z ) 



and is associated with the L-operator 



n-l 

E 

fe=0 



) ® ek,k 



iw z [ei,i (g) ek-i,k + 62,2 ® ek+i,k 



]}• 



This operator is the same as that presented in [18] and is a special case of the i-operator presented 
in [8] , appearing in connection to the chiral Potts model. The former of these articles discusses the 
alteration of the underlying algebraic structure caused by taking the required limit. 
These operators satisfy the Yang-Baxter like relations 



ri2{x)Li3{xy)L23iy) 
Iini2;ri2(2:) Li3{x)L*23{y) 

-Z^i2 )Li3 (xii/i)i?23 (2/1 , 2/2) 
Ll2{x2)Ll3{x2y2)R23{yi, 2/2) 



L23iy)Li3ixy)ri2{x), 
L2z{y)Li:i{x) limzri2(2) 
R2z{yi,y2)Li3{xiyi)Li2{xi), 

R23{yi,y2)L*13{X2y2)Ll2{X2)- 



These relations can be used to construct commuting transfer matrices along with fusion relations. 

In [10] fusion relations were presented for the original DiD^) limit of the Fateev-Zamolodchikov 
model. Following their work we define the vectors 

1 



Vi. = 



V2 

1 

V2 



w 



2 Ul Vk-l ± w 



V2 ffe-1 



~V2 (8) Vk 



and 



Each set of vectors, and {uf}'^^■^, are orthonormal and satisfy 

L{z)v'^ = —iw~^[z ± i\u^ which implies U~^L{z)V = —iw~^ ^ 



[z + z]/„ 




i]In 



10 



for certain matrices U and V where /„ is the identity n by n matrix. 

Before we derive the fusion relations it useful to rescale the i?-matrix, yielding 



R{z)=N{zi,Z2) J2 

a,i,j=l 



L6=l 

where 

n-l 

Nizu ^2) = - n (^1 - ^^')(^2 - w-^'^). 



n 

k=l 



This has the advantage that the entries in R{z) are now polynomial. This rescaling does not effect 
the global Hamiltonian due to the requirement that its trace is zero. 

Periodic boundary models 

Looking at the Yang-Baxter like relations we find Li3{izi)R23{zi, Z2) can be made lower block 
triangular. Algebraically we have found 

V,2 L,,{tZ,)R2s{z,,Z2)V,2 = ^ _ 1)^(^2,^^,^) 

where * represents an unknown operator. This relation can be used to define a functional relation 
for models with periodic boundary conditions. 

t^'^\izi)t^'^Hzu Z2) = (w-^zi + lftS^'>{w-^zi,Z2) + {zi~lft^^'>{w^zi,Z2), where 
t^^^z) = tro[Lo£W-ioiW], 
t^^\z) = tvo[Roc{z)...Roi{z)]. 

Here we have that the functional coefficients of the transfer matrices are only dependant upon z\. 
It is possible to construct a second function relation where the functional coefficients are dependant 
upon Z2- In this case this is achieved by taking the complex conjugate of the above relation while 
noting that 

't(^\z^,Z2)]* =t^''\z*2,zl). 

For more general t'^-'(z), i.e. with twisted boundaries, this will not be case and the two functional 
relations which can be constructed will not be trivially related. Instead of using another 

transfer matrix would need to be constructed using L*{z). 

Using the commuting nature of the transfer matrices wc; write the functional relation as 

A(jzi)A(zi,Z2) = {w^'^zi + l)^A(w"^zi,Z2) + {zi - l)^A(w^2;i,Z2) 

where A(z) and A(z) arc eigenvalues of t^^^ [z) and t^^^ (5) respectively. Following the work of [10, 19] 
we use the ansatz, A(z) is a constant or 

K{z) = c JJ (2:1 - iwyi^h) n (^2 + «w~^y2,fc). (13) 

k=l k=l 

where di, (i2 € N and yi.k, j/2,fe € C. Dividing the functional relation through by A(^i, 2:2) and taking 
the limit zi — )■ iwyij it follows by looking at the residues that 

( 1 + = f\( ^I'fe ~ ^^y^'i 

V 1-iwyij J lj[\yi,k - w-'^yij 

for I < j < di. It is likewise possible to construct a second Bethe equation from the second func- 
tional equation that was mentioned. The second equation in this case can be obtained by replacing 
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yi^k with 2/2, and di with in the above equation. 
Open and braided closed boundary models 

To construct fusion relations for models with open or braided closed boundary conditions we addi- 
tionally require the operator 

n-l 

L'{z) = ^ {[ei,i (g) ek+i,k + 62,2 <S) ek-i,k] - iw~'^z (w^''ei,2 + t«~^''e2,i) ® ek,k] oc 
as L{z) does not satisfy unitarity. It is again possible to diagonalise this operator using the property, 
L\z)u^ =Ti[z±i]vt which implies y-^L'(z)t/ = z f + '1^" V 

Thus is possible to diagonalise both L{z) and L'{z) such that corresponding blocks have the same 
zeros. This allows the construction of fusion relations for both open boundary and braided closed 
boundary conditions. 

One example of functional relations for the open model is 

t(2)(izi)t(3)(zi,z2) = f{w-^zx)t'^^\w-^zuZ2)+g{w-^zi)t^^\w''zi,Z2), whcrc 
= tvo[Loc{z)...Loi{z)L'o^{z)...L'^c{^)], 
t^^\z) = tvo[RociS)-RoiiS)Rio{S)...Rci{z)] 



g{z) = (-1) _ [wz-l\ , 

Here t^^^ (z) will produce an open quantum chain with non-interacting boundary terms. Again taking 
the conjugate equation of this equation yields and additional functional relation whose functional 
coefficients depend upon Z2. Incorporating interacting boundary terms will make the two func- 
tional relations non-trivially related. Using the ansatz given in Equation (13) for the eigenvalues of 
t^^\zi,Z2) yields the Bethe equations 



( 1 - ^'yjj ] ( ^ + w-^vIj \ f l + iw-^yg^X " ^ ^ / ya,k - 



w '^ya,j 

for a G {1, 2} and 1 < j < da- 

We also consider the braided closed boundary models. It is possible to consider the four different 
possibilities outline previously. We let 

zo = (^01, ^02) e {(0, 0), (0, 00), (00, 0), (00, 00)} 

and define the operators 

L(z) = lim ( —^L'(x)] and .R = lim (l + zf^) ^(l + z7^) ^ R2i(z) 

We also make use of the following constants 

60 = 1, ^'00 = —iw~^, Co = 1 and Coo = i'w~^. 
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With these constants we have the following functional relations, 

t^''''\iZi)t^^\zuZ2) = bi^{w-^Zi + lft^''\w-^ZuZ2) + ci^izi-lft^^\w^Zi,Z2), 
t^^'^\-iZ2)t^''\zuZ2) = {KJ''{W\2 + l)''t^^\zWz2) + (0^(^2 - l)^*^^^'^! , «^"'^2), 

^(2,i)(-^^ = tiQ[LQc{z)...LQi{z)Loi{zoi)...LQc{zQi)\ , 

t'-^'^Hz) = tro [L;^^(z)...L^i(z)ISi(^02)...LS£(^02)] , 

= tro [Eo/;(5)...i?oi(^)^oi-^i£] • 

Here we note that the two functional equations arc non-trivially related. Using these functional 
relations and the ansatz presented in equation (13) we have the Bethe equation 

\CzoaJ V ^-iwyaj J l}[\ya,k - W-^ya,j J ' 

for a e {1,2} and 1 < j < da- These equations at most differ by a factor of (—1)^ on the left hand 
side of the equation compared to the periodic boundary case. 



5 Summary 

Taking a special case of the Fateev-Zamolodchikov model yields an i?-matrix which intertwines the 
coproduct of D{Dn). Subsequently models with periodic, open or braided closed boundary condi- 
tions will have the underlying symmetry of D{Dn) (or one of its subalgebras). The Hamiltonians 
presented are seen to be composed of two Hamiltonians with a free coupling parameter. Functional 
equations were also constructed through fusion which allowed for Betlic equations to be determined. 
In the simplest case of n = 3 the ground state energies of the periodic and open boundary models 
have been explicitly computed [19] along with one case of the braided closed boundary models. 
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